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Using two flavors of maximally twisted mass fermions, we calculate the S-wave pion-pion scat- 
tering length in the isospin 1 = 2 channel and the P-wave pion-pion scattering phase in the isospin 
7=1 channel. In the former channel, the lattice calculations are performed at pion masses rang- 
ing from 270 MeV to 485 MeV. We use chiral perturbation theory at next-to-leading order to 
extrapolate our results. At the physical pion mass, we find m^a 1 ^ = —0.04385 (28) (38) for the 
scattering length. In the latter channel, the calculation is currently performed at a single pion 
mass of 391 MeV. Making use of finite size methods, we evaluate the scattering phase in both the 
center of mass frame and the moving frame. The effective range formula is employed to fit our 
results, from which the rho resonance mass and decay width are evaluated. 
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1. Introduction 

Experimentally, hadron-hadron scattering is an important method to study the strong inter- 
actions. Among the various scattering possibilities, pion-pion scattering is the simplest and best 
understood one due to the fact that the underlying chiral symmetry strongly determines the low 
energy behavior of the pion-pion scattering amplitude. Despite its simplicity, pion-pion scattering 
offers us a lot of information on the strong interactions. In the isospin 1 = 2 channel, near threshold 
the S-wave scattering length is used to determine the corresponding low energy constants (LECs) 
of chiral perturbation theory (#PT). In the / = 1 channel, the prominence of the rho resonance is 
clearly observed. By measuring the P-wave scattering phase, the parameters for the resonance mass 
and decay width can be extracted. In the 1 = channel, the sigma resonance appears in pion-pion 
scattering. In contrast to the rho resonance, a precise identification of the sigma resonance remains 
a great challenge because the large decay width of the sigma causes a strong overlap between it 
and its background. 

Pion-pion scattering is non-perturbative in nature at low energies. Therefore, it should be stud- 
ied with a non-perturbative method like lattice QCD. In the center-of-mass frame (CMF), a direct 
lattice QCD determination of the scattering phase is possible by employing Liischer's finite-size 
methods [|], |[ |||], which establish relations between the discrete energy spectrum in the finite 
volume and the elastic scattering phase in the infinite volume. In the moving frame (MF), where the 
total momentum of the pion-pion scattering system is fixed to be a non-zero value, one can evaluate 
the scattering phase by using the method proposed by Rummukainen and Gottlieb [g], which is an 
extension of Liischer's method to MFs. To perform our calculations, we use the Nf = 2 maximally 
twisted mass fermion ensembles from the European Twisted Mass Collaboration (ETMC). Due to 
the properties of twisted mass fermions at maximal twist, our calculation is automatically accurate 
to 0(a 2 ) in the lattice spacing, a. 

In this paper, we present a calculation of the S-wave pion-pion scattering length in the 1 = 2 
channel and the P-wave scattering phase in the 1 = 1 channel. A calculation of pion-pion scattering 
in the 1 = channel using 2+1 flavors of domain wall fermions has been reported recently [|7p. Al- 
though the object of our investigation is simply pion-pion scattering, the approach to study scatter- 
ing from finite size methods in lattice QCD is universal and can be applied to other meson-meson, 
meson-baryon and baryon-baryon scattering systems. 



2. 1 = 2 channel 



In the 1 = 2 channel, the lattice calculation is performed in the CMF. As mentioned in the 
introduction, Liischer's finite size method relates the energy levels of two pion states in a finite 
volume to the scattering phase in the infinite volume. For the case of two pions with zero relative 
three-momentum, this method establishes a relationship between the lowest energy eigenvalue E 
in a finite box of size L and the corresponding scattering length a 1 ^ 2 [0] : 
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where c 



-2.837297 and C2 = 6.375183 are numerical constants. Thus the above result allows 



us to convert a lattice determination of the energy shift dE 1 ^ 2 into a calculation of a 1 ^. 
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To extract the energy shift, SE 1 ^ 2 , and then the scattering length, a 1 ^ 2 , we use the two flavor 
maximally twisted mass fermion configurations from ETMC [[§, |9[ FOfl. 1 The pion masses range 
from m % = 270 MeV to 485 MeV. For most of the ensembles, the lattice spacing is a = 0.086 fm 
and the box size is L = 2.1 fm. For the lower pion masses the volume is increased to L = 2.7 fm. 
Additionally, we perform a check for lattice artifacts with a single calculation at a finer lattice 
spacing of a = 0.067 fm. All the results for the scattering length are shown in Fig. [l]. 
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Figure 1: Chiral extrapolation for the 1=2 pion- 
pion scattering length. The results in this work 
are shown together with the lattice calculations of 
NPLQCD [O 111 and CP-PACS rfl4|| and the ex- 



perimental data from E865 at BNL 
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Figure 2: Difference between the lattice calculation 
of the scattering lengths and the LO ^PT prediction. 
The scattering lengths agree statistically with the LO 
jPT prediction for pion masses ranging from m K = 
270 MeV to 485 MeV. 



The next step is to extrapolate the scattering length to the physical limit. Here, we make use of 
NLO #PT for the pion-pion scattering length, which has recently been studied in the twisted mass 
case [16]. The #PT fit curves are shown in Fig. [l| In the same figure, we also provide a compari- 
son to the lattice results of NPLQCD [fil, ph and CP-PACS [PI] and the experimental data from 



E865 at BNL [15]. We find general agreement between our calculation and the results of NPLQCD 
at similar pion masses. Additionally, we find agreement with the experimental determination of 



m n a^ n . To highlight the impact of the NLO terms in the %YI description of the pion mass depen- 
dence of rnxolj^jr, we show the difference between the lattice calculations of the scattering length 
and the LO ^PT prediction in Fig. ||. We find that the scattering lengths statistically agree with 
the LO #PT result for all lattice calculations with m n < 500 MeV. Accordingly, the NLO %PT 
functional form provides a reasonable description of the lattice results in the same region of m n . At 
the physical pion mass, we obtain the final results 



0.04385 (28) (38) and l l ~ 2 {ii = f KjPby ) = 4.65 (.85)(1.07) , 



where the first error is statistical and the second is our estimate of several systematic effects. For 
more details, we refer the reader to our recent paper [|lt|]. 



The role of the neutral pion in our setup is discussed in some detail in Ref. [1 10 
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3. / = 1 channel 

In the 1=1 channel, the rho meson decays into two pions in the P-wave. As the case of S- 
wave, one can make use of finite size methods to calculate the P-wave scattering phase. However, 
the extraction of the rho resonance from the scattering phase is non-trivial for several reasons. First, 
only when the pion masses are light enough to satisfy the requirement of m n < m p /2, is it possible 
for the rho to decay into two pions. Second, the standard form of Liischer's method is derived to 
address the elastic scattering process, so the interesting energy spectrum should be smaller than 
\m n to avoid the inelastic scattering. Third, because of the finite volume, the energy spectrum of 
pion-pion scattering states is discrete, which translates into scattering phases at discrete energies. 
Therefore, an analytic expression of the scattering phase is required to describe its dependence on 
the energy spectrum. Usually, one employs the effective range formula to meet this demand: 



,2 



I? 



tan 8i (k) = ^ — , =— r- , k = \ Er M /A — mL 
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where 8i (k) is a P-wave scattering phase in the 1=1 channel and Eqm is the center-of-mass energy. 
In the MF, Eqm is simply given by E CM = E 2 — P 2 , where E is the discrete energy eigenvalue and 
P is the total momentum of the MF. 2 Thus, in the effective range formula, only two parameters 
are undetermined, Mr and g p%K , where Mr denotes the resonance mass and g pn]t is the effective 
p — > 7171 coupling constant, which largely determines the size of resonance decay width: 

„ Snnit k^ 



6n M 2 ' 
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By fitting the discrete scattering phases to the effective range formula, one can evaluate the pa- 



rameters Mr and g pK7[ and then determine Fr. Conversely, by using the latest PDG [ |18| ] values of 
«% = 139.5702(4) MeV, M p = 775.49(34) MeV and F p = 149.1(8) MeV, one can also evaluate 
gpnn = 5.98(2) at the physical pion mass. 

With the effort required to simulate with light up and down quark masses, the condition of 
m % < m p /2 has been satisfied by only a few lattice calculations [^, |2(], 21]. So far, all these 



studies concentrated on one or two scattering phases for each ensemble. Since the effective range 
formula carries two unknown parameters, more scattering phases are needed for a precise fit. To 
accomplish this goal, a natural way is to calculate the energy spectrum of the higher excited states. 
In order to isolate the ground state and the first excited state, we set up a 2 x 2 correlation matrix 

/ ((^)t( ? )(^)(0)) <(«*)t(f)p(0)) \ 

C2x2W ^ <pt(0(**)(o)) <pt(Op(o)> ) ' {5A) 

where the interpolating operator (7i7i)(t) has the same quantum numbers, J PC = 1 , as the inter- 



polating operator p(t). After diagonalization of the matrix in 3.1, we obtain the energy eigenvalues 
of the ground state and the first excited state, E n (n = 1 , 2), and then convert them into the scattering 
phases. More ambitiously, constructing a. N x N matrix allows us to look at even higher excited 
states. However, the realistic computation of the scattering phase at higher energy, E n (n > 2), 

2 To reduce lattice discretization effects, we use the relation cosh(£'c^) = cosh(£') — 2sin 2 (P/2) instead. 
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remains a challenge due to the poor signal-to-noise ratio and the restriction of E" < 4m n . Another 
way to determine the scattering phase at more energies is to perform a lattice calculation in the 
MF. In our case, we use a MF with a total momentum P = e^ln/L, which provides us another two 
scattering phase points. In principle, by performing the lattice calculation in the MF with other 
total momenta, for example P = (e\ + e2)2n/L, it is possible for us to collect even more points. 
However, we must be careful in choosing the MF because in some MFs the ground state and the 
first excited state are nearly degenerate and isolating them becomes very difficult. 
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Figure 3: Center-of-mass energies in the CMF and 
MF. The energies evaluated from the correlation ma- 
trix are compared with the ones evaluated from the 
diagonal mattix elements. 



Figure 4: Four scattering phases calculated on the 
lattice together with the effective range formula fit. 
At the position where the scattering phase passes 
Tt/2, the resonance mass Mr is determined. 



In this work, both the calculations in the CMF and MF are performed by using Nf = 2 dynam- 
ical maximally twisted mass fermions. The corresponding lattice parameters are m % = 391 MeV, 
a = 0.086 fm, L = 2.1 fm and m n /m p = 0.4. By diagonalizing the 2x2 matrix in [3.1| , the Eqm 
of the ground state and the first excited state are evaluated and shown in Fig. [| In order to inves- 
tigate the effect of diagonalization, we also perform a study that utilizes only the diagonal matrix 
elements. We see that in the CMF, there is a strong mixing between the ground state and the first 
excited state in ((nnY (t)(nn)(0)) . While in the MF, a similar situation happens to (p^(t)p(O)). 
So no operator safely provides us the ground state energy in both frames. Therefore, introducing 
the diagonalization method to the calculation of rho decay becomes essential. 

As shown in Fig. ||, all the four values of Eqm are smaller than 4m K . Unlike continuum QCD, 
twisted mass LQCD violates the symmetries of isospin and parity. As a result, it is possible for 
the rho to decay into three pions, which means that at non-zero lattice spacing the upper bound 
of the elastic scattering region is lowered to 3m n . Additionally, the isospin symmetry breaking 
causes a mixing between the 1=1 channel and the possible 1 = and 1 = 2 channels. Thus, a 
corresponding modification would be required to adapt Liischer's method to the isospin mixing 
case. In Ref. |0], a significant effort was made by us to attempt to find these effects in the 1 = 2 
channel, but no compelling evidence was found. However, the effects of isospin violation in the 1=1 



channel are expected to be more severe [16]. Just for the purposes of these proceedings, we assume 
that such effects are small. Therefore, we convert all the four Eqm values into the corresponding 



scattering phases using the normal method. As in Ref. [ 17], we will eventually examine the effects 
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Figure 5: L dependence of E£ M (n=l,2) at a pion 
mass of m, = 391 MeV, using values of Mr and g pnn 



in Eq. 3.2 



Figure 6: L dependence of E£ M (n=l,2) in the phys- 
ical limit, defined by the latest PDG [ |T8| | values of 
Mr and g pmc . 



of parity breaking carefully and complement the current calculation with a calculation at a finer 
lattice spacing to explicitly check for any strong lattice artifacts for / = 1. 

The results for the scattering phase are shown in Fig. ^| together with the effective range for- 
mula fit. At the position where the scattering phase passes % /2, the resonance mass Mr is deter- 
mined. Additionally, the values of g p7[K and Fr are also evaluated from the fit. Our final results 
are 



aM R = 0.4186(56) , g p%% = 6.16(48) and aT R = 0.0217(44) . 



(3.2) 



Here, our result for g p7t7[ at m n = 391 MeV agrees statistically with that at the physical pion mass, 
which hints that the pion mass dependence of g p7[K might be weak. However, we can not make 
any strong statements here since our calculation of rho decay is only performed at one pion mass 
and our errors are rather large. To determine the chiral limit of g pnn and Mr, we will perform our 
calculation at smaller pion masses. 

Usually one takes the ground state energy E^ M in the CMF as the rho mass. This is correct 
in the case m n >m p jl where the rho is still a stable particle and E^ M has only an exponentially 
suppressed L dependence. However, when m n becomes smaller the rho becomes unstable and E^ M 
begins to depend more strongly on L. As an example, we calculate Ecm{L) by combining Luscher's 



formula and the effective range formula with the parameters Mr and g p7t7l given in Eq. \i.2[ As 
shown in Fig. ||, the L dependence of the lowest level is visible but still weak at a pion mass of 
m K = 391 MeV. Decreasing m n further, more phase space becomes available for the rho to decay 
into two pions. Thus, assuming g pK7[ is roughly constant, the decay width becomes larger and the 
lowest level in the CMF begins to look more like a scattering state. Fig. ^ shows Ecm{L) in the 
extreme limit at the physical pion mass. There we see that E^ M drops so rapidly with L that it bears 
no resemblance to a stable state with mass Mr. 



4. Conclusions 

We have calculated the S-wave pion-pion scattering length in the 1 = 2 channel and the P-wave 
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pion-pion scattering phase in the 7 = 1 channel using Nf = 2 maximally twisted mass fermions. In 
the former channel, the pion masses ranged from 270 MeV to 485 MeV. Using #PT at NLO, we 
extrapolated our results for the scattering length to the physical limit, where we found m^a 1 ^ = 
-0.04385 (28)(38) and Z^= 2 (jU = f„ ;phy ) = 4.65 (.85)(1.07). In the latter channel, we performed a 
calculation at m n = 391 MeV, a = 0.086 fm and m n /m p = 0.4. Making use of finite size methods, 
we evaluated the scattering phase at four energies, to which we fit the effective range formula and 
found the results aM R = 0.4186(56), g pmt = 6.16(48) and aT R = 0.0217(44). 
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